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!>• ' Abstract 

o\ , 

, This paper deals with the comparison between the strong Thomas- 

Fermi theory and the quantum mechanical ground state energy of a large 
(3 ' atom confined to lowest Landau band wave functions. Using the tools of 

microlocal semiclassical spectral asymptotics we derive precise error esti- 
mates. The approach presented in this paper suggests the definition of 
a modified strong Thomas-Fermi functional, where the main modification 
consists in replacing the integration over the variables perpendicular to the 
magnetic field by an expansion in angular momentum eigenfunctions. The 
^ ' resulting DSTF theory is studied in detail in the second part of the paper. 

> 

■ 1 Introduction 

H , 

5^ , In this paper we study semiclassical theories describing the ground state energies 

of heavy atoms in strong homogeneous magnetic fields, where additionally the 
electrons are confined to the lowest Landau band. 

An atom with N electrons of charge — e and mass rrie and nuclear charge Ze 
is described by the nonrelativistic Pauli Hamiltonian operator 



^<j<N ^ \ j\ J l<i<j<N ' 



acting on the Hilbertspace L^(M'^, C^) of electrons. The units are chosen 

such that h = 2me = e = 1. The magnetic field is B = (0,0,5), with vector 
potential A = ^B{—X2,Xi,0), where B is the magnitude of the field in units of 

2 3 

Bo = "^"^ = 2.35 ■ lO^Gauss, the field strength for which the cyclotron radius 
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l-B = {he/ {eB))^/"^ is equal to the Bohr radius oq = /{rriee^). The ground state 
energy is 

E^{N, Z, B) = inf {(t/^, Hni/j) : G domain H^, ijj) = 1}. (1.2) 

Recall that the spectrum of the free Pauli Hamiltonian on £^(M'^;C^) for one 
electron in the magnetic field B, 

Ha = [(T ■ i-zV + A{^))]\ (1.3) 

is given by 

pI + 2uB z/ = 0,1,2,..., p,eR. (1.4) 

The projector IIo onto the lowest Landau band, z/ = 0, is represented by the 
kernel 

no(x, x') = ^ exp |^(xx X xl) ■ B - ^(x^ - x^^fij 5(2 - ^')A> (1-5) 

where x^ and z are the components of x perpendicular and parallel to the mag- 
netic field, and P| denotes the projection onto the spin-down component. 
In this paper we are especially interested in the ground state energy, 

e2,UN,Z,B) = irif (^,<i7^<^), (1.6) 

||V'||=i 

where 11^ denotes the iV-th tensorial power of IIo. Lieb, Solovej and Yngvason 
pointed out that for B ^ Z^^^ the electrons are to the leading order confined to 
the lowest Landau band, which is expressed by the following theorem. 



1.1. THEOREM. (^SY^, Theorem 1.2) For any fixed X = N/Z there is a 
6{x) with 6{x) — > as X — >• cx) such that 

^clf > > ^If (1 + S{B/Z'/-')) . (1.7) 

The energy ( |1.6[ ) can be approximated by means of the STF-functional (Strong 
Thomas- Fermi) 

^'^"[p] = ^/p^-/ Vp + Dip,p), (1.8) 

V{x) = Z/\x\ and D{p, p) = i(p, \x\-^*p). In |LSY2i it is shown that EQ/E^tf ^ 
1 if Z ^ 00, B/Z^ and B/Z^/'^ 00, where 

E'^^iN, Z, B) = MiS'-'^lp] \p>0,pe D'^^, Ip<N}, (1.9) 

with an appropriately chosen domain D^'^^. Combined with Theorem LI this 
implies Theorem |1.2| . 

1.2. THEOREM. (^SY^) If Z 00 with N/Z fixed, B/Z^ andB/Z^/^ - 
00, then 

Elf (iV, Z, 5)/ESTF(iV, Z, B) ^ 1. (1.10) 
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1.1 Comparing the STF energy with the QM ground state 
energy in the lowest Landau band 



In this paper we not only want to give a direct proof of Theorem p. .21 , but we 
want to derive precise error estimates. In this respect our procedure is related to 
and [|T|. Our main theorem is the following: 



1.3. THEOREM. Let N ^ Z and Z^l^ < B < Z'^ . Then 

\e2,^,{N,Z,B)- E'^^''{N,Z,B)\ < CB^'^'Z''"' (1.11) 
for an appropriate constant C . 

1.4. Remark. In ||n|] Ivrii estimated the difference between the full quantum 
mechanical energy and the MTF energy, which is given by minimizing the 
MTF functional (cf. [ LSY2|| ), where all Landau levels are taken into account. 



The main theorem in |^ reads: 

1.5. THEOREM, (g^ Theorem 0.2) Let B < Z^ and N ^ Z, then 

\E'^{N, Z, B) - E^^'^{N, Z, B) -^Z^\<Ri + R2, (1.12) 

with 

Ri = CZ^'^{N + Bf'^ and R2 = CZ'^I^B^'^. (1.13) 

1.6. Remark. Although true for all B < Z^, it should be noted that only for 
B < Z^^ are the error terms Ri and R2 smaller than the Scott term jZ . 



One of the main difficulties Ivrii has to cope with in the proof of Theorem |1.5| is 
the fact that the self-consistent MTF potential is not smooth, because it includes 
all Landau levels. So he has to create an approximating C°° potential in order 
to apply the tools of microlocal semiclassical spectral asymptotics. Fortunately, 
in our case of Theorem |1.3| we need not care about such problems, since the 
STF potential (see (|1.14|) ) has all required properties for semiclassical spectral 
asymptotics. 



Moreover, in Theorem |1.5| Ivrii already captures ( |1 . 1 1| ) on the region where 



only the lowest Landau band is occupied, i.e on {x||x| > CqZ/B} with a large 
constant Cq. 



We mention some important steps of the proof of Theorem \L3 . 
Let (j)^"^^ denote the effective STF potential (for simplicity think of the neutral 
case N = Z) 



ST^^-^-Z|xr^-pS^%|x|-\ (1.14) 



[X 
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where p^"^^ is the minimizer of the STF functional ( |1.8| ). We will see in Section 
2.1 that the main contribution of ( p. . 1 Ij ) is given by 



Tr[no(ifA-0'^''(x))no]- 



B 
2^ 



;i.i5) 



with = min{0,)f:}. Recall 



E 



STF 



B 
2^ 



Zn 



1.16) 



With the decomposition C^{M.^, rfx; C^) = C^{M.'^, rfx^) ® C^iR, dz) ® the pro- 
jector Ho can be written as 



[in) 



m>0 



where 0^ denotes the function in the lowest Landau band with angular momen- 
tum —m < 0, i.e., using polar coordinates {r,ip), 



B 1 fBr 



2\ rn/2 



27r ym,! 

Using this and Hx^m = 0, we can write 



;i.i8) 



;i.i9) 



m>0 



By means of the above decompositions one gets the relation (cf. 0] Theorem 
3.13) 



with 



0^^= / rfx^0"'^'^'(x)|0™(x^ 



STF/ 



;i.2i) 



Next we multiply the m-th term of the right hand side of ( |1.20| ) with ^Xm(x^ 
where 

1 for ^2m/B < |x^| < ^2{m + l)/B 
otherwise, 



1.22) 



and integrate over x^, which is just an identity operation. Since we are allowed 
to put the sum into the trace as well into the []_ bracket we arrive at 

Ti[M-dl - 0'^^(x)]no]- = d^^TTcH^)[-dl - 0^(x)]_, (1.23) 
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with 



kSTF 



;x) = 5^Xn.(xx)0rW- 



1.24) 



The Equation ( |1.23| ) follows from the fact that the terms Xm{'^±] 
have disjoint supports. 

Hence, ( p..l5| ) can be written as 



B 
2^ 



R2 



dzdp 
27r 



,STF, 



x)]- 



We shall estimate ( |1.25| ) by splitting into the following two terms: 



^ / dx^(Tr^.(M)[-9,^-0STF(x)]_ 



dzdp 

271 



kSTFi 



x)]- 



dr^) 



1.25) 



1.26) 



and 



2tT /in.2 



^[^[p2_0STF(x)]__ 

2tt 



(1.27) 



1.2 Modified STF functionals 



From Equations (|1.20 ) and ( 1.23|) it is apparent that the STF energy (|1.9| ) is not 
the most natural semi classical approximation of -E^^f. As already argued in 



(|1.20|) suggests the definition of a functional, where the integration over x^, the 



variables orthogonal to the magnetic field, is replaced by an expansion in angular 
momentum eigenf unctions in the lowest Landau band. This leads to a discrete 
STF functional (DSTF) depending on a sequence of one-dimensional densities 
P = (Pn)neNo, i-e. 

^DSTF[^] =J2U[ Pm{zf-Z I V^{z)p^{z)dz\ +D{p,p), (1.28) 

where k = 7r^/3, 

-D(P, P) = ^^ / Vrn,n{z - Z)prn{z)pn{z')dzdz' , (1.29) 
m,n 

and the potentials Vm and Vm,n are given by 

Vm{z) = J j^|0m(x_L)p(ix_L, 

^„,„(,_,') . / l^"'"^)l'l^f' ^'I' dx.dxl. (1.30) 
•/ |x X I 
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An equivalent functional, depending on a three-dimensional density p, can be 
obtained as in (|1.23|) , if in STF theory the Coulomb potential is replaced by 



|x| ^ = V'Xm(x_L) / rfx_L|x| Vm(x_L). (1.31) 

The resulting modified STF functional is given by 

^""'^"[P] = ^ / ^^'^W' - / ^xW^V(x) + D{p,p), (1.32) 

Dip,p) = ^^ J dxdyVm,n{z - z')xmi^±)Xn{y±)pi^)piy)- (1-33) 



with 



Let p^^^'^ = {p^^^^{z))m and pi^s^F ^ pMSTF^^^^ ^.j^g ground state densities of 
(|1 . 28|) respectively ( |1.32|) corresponding to a fixed particle number A^. Then the 
relationship between the densities reads 

P^^^'^(x) = ^Ep!^'''(-)^-(-^)- (1-34) 

m 

Furthermore the energies are equal, 

^MSTF^^^ Z, B) = E°STF(Ar, z, B). (1.35) 

Since a short computation shows that for B > Z^^^ the difference between the 
D(M)STF energy and the STF energy is smaller than B'^^^Z^^^ the estimate ( |1 . 1 1| ) 
with STF replaced by D(M)STF immediately follows for this region (B > Z^^^). 

1.7. THEOREM. Let Z^/^ < B < Z^ and N ^ Z . Then 

\e2,^,{N,Z,B) - E''^^^'^^^iN,Z,B)\ < CB^"'Z^I\ (1.36) 

2 Proof of Theorem |l73l 

2.1 Derivation of lower and upper bounds to 



First of all, recall that the TF equation satisfied by the minimizer of (|1.8|) under 
the constraint j p = N is 



K{p''y/B'' = [Z|x|-i - p^''' + = + (2.1) 

where u = v{N) is the chemical potential corresponding to the electron number 
A^. Using (|2.1|) one sees that the STF energy (|1.9|) can be written as 



E^^^ = ^ [ d^J - 0^T^(x) - + uN- Dip'^^,p'^^). (2.2) 

27r ./iB2 Its Zn 
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This expression will be used for deriving upper and lower bounds to the quantum 
mechanical ground state energy -E'^^f- 
Lower bound: 



Let ip denote a ground state wave function of if at, then we can write -E^^f 



E. 



N 



i=l 

+Ni^ - 2D{p^, pSTF) ^ |x, - x,|-V), (2.3) 

i<j 

where we have added and subtracted the term p^"^^ * |x|^^ — u and used the 
definition 



By means of the Lieb-Oxford inequality ||L0| 



El-/-. 



"V)>C((J#.P#)-1.68 / p*' 



i<j 



4/3 



(2.4) 



(2.5) 



^ can be bounded from below by 

^If > Tr[no(i/A - 0^^" - ^)no]- - D{p^^^, pSTF) ^ _ i_Q8 I ^4/3^ ^2.6) 



where we have used that 

D(p,-pS^^p,-p^^^)>0. 
Furthermore by (|1.23|) and ( p.2|) we get 



STF 



7^l - 1.68 J pf, 



(2.7) 



(2i 



with 



7^l 



2^ 



(ix_L (^Tr£2(iR)[-a 



kSTF( 



iSTF/ 



(2.9) 



Since is a. ground state wave function, or at least an approximate ground state 
wave function, we can estimate (cf. |[LSY1|| (8.5)) 

pf < ( j pIY'\ I p,r' < const.{B^\E^^^\Y''N'/' 

< constZ'/'N'^/^'B^/' < CB^"'Z^'\ (2.10) 
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using N r-. Z and B> Z^/^. 
Upper bound: 

For every fixed integer N and a normalized iV-particle wave function ip we liave 



N 



i=l 



+Nu - 2D{p^, pSTF) ^ |x, - x^rV). (2.11) 

We set 

^ = -^0iA...A0Ar, (2.12) 

wliere (pi is tlie eigenvector corresponding to the i-th lowest eigenvalue Aj of the 
one-particle operator 

IIo{Ha - Z\^\-' + pST^ * |x|-i)no. (2.13) 
By means of the decomposition 

N 

Tr[no(i/A(x)-0S^^(x)-z/)no]- = 5^[A.-z/]_ = 5^(A.-z/)+ (A.-z/), 

i 1=1 \fi<\i<v 

(2.14) 

and the equation ( |2.7|) we can estimate (|2.11|) as 

^Snf < Tr[no(//^-0^^^-z.)no]_ + z/iV-D(pST^pST^) 

+D(p^-pST^p^-pS^^)- 5^ (A,-z.), (2.15) 



Ajv<Ai<!/ 



which implies 
with 



^clf < ^1 + ^2 + i?3, (2.16) 

i?2 =- J] (A,-z/), (2.17) 

Ajv<Ai<i/ 

i?3 =I^(p^-p'^^p^-p'^^). (2.18) 

Since it is difficult to tackle directly the term i?2, we estimate R2 by \\n — u\ 
multiphed by the number of eigenvalues of the operator ( ^.13| ) between Aat and 
V. We know 

/" f dzdp ^ ^ 2 jSTF/ \ \ _ f STF 



^ / c^xW =e_(j>^-0^^Mx)-z.)= / p^''=N, (2.19) 
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with 

^^^^^ I otherwise. (2.20) 
So the number of eigenvalues of ( |2.13| ) between Aat and u can be expressed by 

Tre_(no(-9,^-0ST^(x)-z.)no)-^ J d^^j ^e^p'-^P'^'i^)-^)- (2.21) 

Mimicking the derivation of ( |1.23| ), with being replaced by 0-(t), leads to 

(2.21) = ^ £ dx^ (^Tr^.(M)e_(iJ.^ '"^'1 '^^-^^-^ - '^)) ' (2.22) 

where we have defined 

H^^ = -dl - 5s^5(x) and = - 0^^''(x). (2.23) 
Hence, instead of R2 we estimate the error term 

7^2 = IAat - f (ix^ (tic2(r)Q-{H^^ - i^) - f ^^^0-(/ix^ - K 



(2.24) 



We introduce the notation 

e{H^^,fi) = Q_iH^^-fi), (2.25) 

the projector of the operator i^xx onto the eigenspace corresponding to the eigen- 
values smaller or equal to /i. Let ip be given by ( |2.12| ), then we have 

B ~ 

P^(x) = — e(z,2;;ifxx,-^7v), (2-26) 
zvr 

whereas the semiclassical density (F^^ can be written as 
Furthermore we introduce the auxiliary density 

p(x) = ^[0STF(x) + A^]V2_ (2.28) 

In order to bound the error term i?3 from above, we combine the two terms 

7^3 =D{p^-p,p^--p), (2.29) 
7^4 = D{p^^^ ~ -p, p^^"" - -pl (2.30) 

which are easier to handle than alone. Observe that by convexity one has 
-R3 < 27^3 + 27^.4. In the next sections we will separately have to carry out the 
estimations 

TZi < CB^I^Z^I^ 'ii = 1, ..,4. (2.31) 
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2.2 Methods used in the proof 

The methods used here in order to estimate the error terms TZi-4, have been 
estabhshed in ||IS| and Q. For sake of better understanding we will state here 
the most important theorem, which we use throughout this section. 
Consider the Schrodinger operator 



H 



Its symbol is denoted by 



2 



on 



(2.32) 



(2.33) 



In addition there are the following conditions imposed on the real potential 
There are Lipschitz functions Z(x) > and /(x) > 0, such that: 



(n) 
(iii) 



|V/(x)| < M, 

c/(y) < /(x) < C/(y), 

|9Xx)| < a/(x)'/(x)-i'^i e N'^ 



(2.34) 
(2.35) 
(2.36) 



Under these assumptions Ivrii and Sigal have proved the following theorem: 
2.1. THEOREM. (/^ Theorem 7.1.) 

Assume conditions (i) — (iii) are obeyed and let ip be smooth and obey \d'^ilj{x.)\ < 
CJ{'x)~^'^^ for any u. Let (?s(A) = [— A]^ for some s E [0, 1]. Then 



dxdp 



27T 



ipgsih) 



< C 



(ixmax 



supptp 



/(x)/(x) 



a—s—d 



(2.37) 
/(x)-2--^ (2.38) 

(2.39) 



The most important tools for the proof of Theorem |2Tl| are multiscale analysis 
and semiclassical spectral asymptotics. First of all, the domain, i.e. the support 
of ipy is covered by a countable number of balls. Then on each of these balls 
i?(y, /(y)), the operator H is transformed into 



Here a = 1 if either d > 2 or d = 1 and (f) obeys 

|0(x)|+/(x)|V0(x)|>5/(x)2 
on {x |/(x)/(x) > 1}, with some £ > 0, and a = 1/2 otherwise. 



K, = f^'umUil)-' = -yA - V(x) 



(2.40) 



with h = l{y) ^/(y) ^, by means of a unitary scaling transformation U{1), which 
maps the ball i?(y,/(y)) into 5(0,1). Next Theorem H 6.1, which we state 
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below, is applied to the sum of the negative eigenvalues of gs{Kh). After rescahng 



and summing over all balls one arrives at Theorem ^TT 
The symbol of Kh is given by 

fc(x,p) =pV2-K(x). (2.41) 

If furthermore all derivatives of V are bounded by a constant, i.e. 

|9^1/(x)|<a on 5(0, 2) Vz, (2.42) 

then the following theorem is valid: 

2.2. THEOREM. Let G C^{B{0,1)) and let be a regular value of the 
function ko restricted to suppip ^ Then for h < 1 

TriijgsiK,)) = R-' j ^^g,{k{x,p)) + 0{h^^'-'). (2.43) 

Assume next that the potential V can be written as 

\/(x) = \/o(x) + /i\/(x,/i), (2.44) 

such that the principal symbol of reads ko{'x,p) = — Vo(x). Then, if 
additionally V fulfills ( |2.42|) (uniformly in h), Theorem p.2| remains valid for 
s = with k replaced by the principal symbol k^. (cf. @ Theorem 4.5.3) 

2.3. THEOREM. Let G C^(5(0,l)) and let be a regular value of the 
function k^ restricted to supptp x R"'. Then for h < 1 

Trii^goiK,)) = h-" [ ^ijg,{ko{x,p)) + Oik'-"). (2.45) 

J 277 

2.3 Relationship between the potentials 0^^^ and 0^^^. 

Next we collect some information about the potentials (p^^^ and (p^^^. 

The scaling functions, which we will use in order to apply Theorem PTTI to the 



operator H^^j^, have to be chosen such that the conditions (p.34|) - (|2.36|) hold for 

the potentials (p^'^^ and (p^'^^, at least away from the origin. 

Since we will see that (p^"^^ behaves like Z|x|~^ for |x| < rs, the edge of the 
STF atom, it is thus natural to define 

/(x) = (const.) |x| and /(x)^ = Z|x|-\ (2.46) 

If we denote the effective STF potential as 

(x) = 0^^^(x) + ^, (2.47) 

then the following Lemma is valid: 
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2.4. LEMMA, (i) (j^SY^) The density p^^^{x) as well as [1/|t^(x)]+ have 
compact support with radius rs < S.Svr 

^1/5^-2/5^ 

(a) For fixed but arbitrary :xi±, (/)^'^^(x) and (f)^'^^ {x.) are G C°°(]R\0) as a function 
of z and 



\d. 



i/ -STF, 



X)| < a/(x)^/(x)-^ |5:0STF(x)| < C:/(X)2/(X) 



(2.48) 



for allu eN and x G . 

Proof, (ii) The property follows from the TF equation ( |2.1| ) and the definition 
([1:241) . Equation (|^ follows from and ([L^ ) and the fact, that 



i5:z|x|-i| < /(x)vr'^— < /(x)2/(x)-^ 



The proof of (i) is given in [LSY2| Theorem 4.11. 



(2.49) 
□ 



2.5. Remark. The estimates ( ^.4(^ ) seem to be very crude, especially in the 
vicinity of rs, but nevertheless they are good enough to provide precise error 
estimates. Next let us try to get an idea how behaves in the vicinity of 

the radius rs- We consider the neutral case N = Z. Since is spherical 

symmetric, we can make the ansatz V^|"'"^(x) = x(r)/r, (|x| = r), which leads by 



O to 

X"{r) = Br^'^x{r) and x(0) = Z. (2.50) 

Around each point tq, this equation has a solution x that can be expanded in 
a series of terms Ci[rQ — r]* with i > 4. In the vicinity of tq = rs we get the 
approximate solution 

Vr-^^[rs-r]l (2.51) 

which shows that tends to as [rs — r]^ as r ^ r^. 

Next we fix a point y G M'^ and we set I = /(y) and / = /(y) = Z^^'^ly]^^'^. 
(We assume that (|2.34|) - (|2.36| ) are fulfilled in i?(y,/(y)). This in our case 
can be done by defining /(y) e.g. as |y|/2.) Furthermore we define the unitary 
transformation 

U{1) : ^(x) ^ /3/2^(^^ ^ ^2.52) 

which maps the ball -B(y, /(y)) to 5(0, 1) and transforms the operator —9^ — 0^"^^ 
into 

-l-^dl-^{lx + y). (2.53) 

Introduce the new potential 

W{x.) = f~^'4F^{lx + y). (2.54) 
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The resulting operator is related to the original one ( |2.23| ) as 

U{l)H^^U{l)-' = fKh, (2.55) 

with 

Kh = -h'^dl - W{-k) and h = {lf)-\ (2.56) 

If we denote W{^) = /~^0^'^^(/x + y), then one easily sees that W can equiva- 
lently be defined by applying ( |1.24| ), i.e. the operation ~ to W, with B replaced 
by B' = Bp. In other words the unitary transformation U{1) scales the magnetic 
field strength B to B' = BP and for the difference W — W we get the following 
Lemma. 

2.6. LEMMA. There exists a function a(x, a) such that 

W{^) - W{^) = aa(x, a), (2.57) 
with a = B~'^/'^l~'^ and a(x, a) fulfills \2.4^ ) uniformly in a for a < 1. 



Proof. Since the potential 0^"""^ is spherical symmetric and Equation ( |2.48| ) is 
fulfilled for derivatives in all directions, we get for x G 5(0, 2) 

\d''W{-K)\ = \f-^d''<f^^{li^ + y)\<C^ for all z/ e (2.58) 

Hence, Vr(x), together with all derivatives, is bounded above by a constant on a 
ball around x = 0. Since the operation ~ smears^ the potential, for every x, over a 
region ~ a in the |x_|_|-direction the difference W — W can be expressed by a times 
a function a(x, a) which is bounded by a constant. Since d'^W{x.) = 9"W^(x) the 
same argument can be given for all derivatives. □ 

Let us rewrite the operator in the form 

K,, = -h'^dl - (H^(x) + aa(x, a)). (2.59) 

In order to be allowed to apply Theorem ^]3| to ( |2.59|) , i.e. in order to guarantee 
that — W{x) is the principal symbol of Kh, it is necessary, that 

a<h ^ B-^/H-^ < l-^f-\ (2.60) 

which leads to the condition 

|y| > Z/B. (2.61) 

Hence, in the sense of Theorem |2.3| , this implies that in the region {x| |x| > 
Z/B} Theorem |2.1| , with s = can be applied to the error terms 7?.2_4, with 
0STF replaced by (p^^^^. Next we decompose into Qi = {x||x| < Z/B} and 
^2 = {x| |x| > Z/B} and estimate the error terms 7^i_4 on each of these regions 
separately. 
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2.4 Analysis in the region Qi 

We first assume that B < Z^. This assumption is made in order to be sure that 
Qi is not completely contained in the non-semiclassical region {x||x| < 1/Z}, 
where each term, the quantum mechanical as well as the semiclassical, has to 
be estimated separately. Furthermore let tp^^^x) be supported in {x| < |x| < 
Z/B{l + e)} and fulfill ^«(x) = 1 in {x| |x| < {Z/B){l-e)}, as well as {d^ip^^^l < 
C„/(x)~" for all neN. 

With respect to TZi and 7^2, we in particular have to estimate the term 

^ I rfx^ (rrc^^^^'^g^iH^^ ~ ""^^ ~ j '^'^'''^ 9s{K^ - , (2.62) 

with ifxx and /ixx given by ( p.23|) . Let h^^ be defined analogously, i.e. hy^^ = 
p2_^j^_ Adding and subtracting / ^iIj^^^ gs{K^- u) in we split p3^ ) 

into 

(2.63) 

and the fully semiclassical part 

R'.m = ^jd^^ (/ '^-^^'''9s{K. ^^^^^^s(v - 

(2.64) 

Since l{x.)f{x.) > 1, which is equivalent to |x| > Z^^, is necessary for being able 
to apply Theorem |2.1| , we have to carry out a corresponding decomposition of 
^li. Let ip^i^ + = ip^^^ be a partition of unity on fii, with 

supp7/>f^ = {x| |x| < Z-\l+e)}, supp^^^^ = {x| Z-\l-e) < |x| < Z/B{l+e)}, 

(2.65) 

and |c?"'?/'j-^'' I < C„/(x)~" for i = 1, 2, with C„ independent of Z and B.. 
2.7. LEMMA. With above definitions we have for ( \2. 6!^) 

Rlii^^^^) < C^i^-^Z'-^^ (2.66) 

Proof. On supp ip^^ we apply Theorem to ( |2.62|) with a = 1 and d = 1, which 
implies for arbitrary but fixed x^ (we may set the chemical potential i> = for 
simplicity, the computations for arbitrary u are essentially the same) 

TTcHM){4'^9siHyJ)- [ ^4'^9sCh.J <C [ rfz/(x)-^-7(x)^ 

(2.67) 
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Hence, multiplying with ^ and integrating over leads to 
2vr J Jz-i 

(2.68) 

In the case of r < the terms of ( p.63|) have to be estimated separately. The 
semiclassical part reads 



271 / 271 



< — [ rfx[0S?5]y2+- ^ CfiZ^^-' (2.69) 
27r Jm<z-^ 



An analogue estimate one derives for ^ J d:x.±Trc2(R){il'\^^ gsiH:)f^j_)) by using [|^] 
Lemma 7.9. □ 



2.8. LEMMA. For ^2M ) we have 



(2.70) 



Proof. Obviously, the main contribution to the magnitude of the semiclassical 
term 



is produced by the Coulomb singularity, i.e 



B 

< — 

27r ./|x|<B-i/2 



(2.71) 



(2.72) 



< CB 



x|<B-i/2 



dx ( — 

r 



s+l/2 



< ^^^+1/2^^/2-1/4^ ^2.73) 



□ 



Hence, we are ready to carry out the estimate of the error terms 7^i_4, re- 
stricted to Qi, which we denote as TZi{Qi). 



2.9. PROPOSITION. For Z^l^ < B < Z^ one has 

n,{n^) = ^ [ dx^ fTr^2(K)(^(i)[#x, - i^]-) - 
< CB^/^'Z'/^ 



271 



(2.74) 



Proof. This is done by putting together Lemma |2.7| and Lemma |2.8| and setting 
s = 1. □ 



Before turning to TZ2{^i) we need a preparing Lemma. 
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2.10. LEMMA. Let be the N-th eigenvalue of ( \2.13l ) and v the chemical 
'potential of belonging to the electron number N. Then 



Proof. We assume now that we have aheady got the estimate 



ASTF^ 



(2.75) 



(2.76) 

which till now we have only proven on Qi, by setting s = in Lemma p.7| and 



Lemma |2.8| . The missing part will be proved in Lemma p.l6| . Since we know by 
definition 



B 
2^ 



B 
2^ 



we get 



^STF _^ ^1-1/2 ) > 



for some u' G [u, X^] . This implies the statement of the lemma. 



(2.77) 

(2.78) 
(2.79) 
□ 



2.11. PROPOSITION. 



7^2(^^l 



Xn - y\ 



B 



27r 
dzdp 
2n 



[ dxxfTr^2(R)(^«e_(i^,^-i/))- 



< CB^/'^Z'/^ 



(2.80) 
(2.81) 



Proof. By Lemma |2.1CI| and combining the estimations of Lemmata |2.7| and ^]8 
with s = 0. □ 



2.12. Remark. We remark here that if one has a partition of unity ipi + (p2 = 
then the relation 



D{f, f) < 2D(M, M) + 2D(M, M) 

is valid, which one gets by the simple inequality 

D{f^l - f^2, f^l - f^2) > 0. 



(2.82) 



(2.83) 



Remark p.l2| justifies the notations 7^3(^2^), 7^4(f2j), since ( |2.82| ) means that 
we can consider each region f2j separately. 



H March 8, 2001 



17 



2.13. PROPOSITION. For Z^l^ < B < we get with 



7^4(^^l) = D(^«(pSTF _ -)^^(i)(^STF _ < CB^/'Z'/'. (2.84) 
Proof. Recall 

P^^^ = f [0^^^ + and -p = [0^^^ + A^]f . (2.85) 
By the Hardy-Littlewood-Sobolev inequality we derive 

7^4(fil) = Dii;^'\p'^^--p),i;^^\p^^^--p))<CB^U^'\p^^^~p)\\l/, 

< CB^lXj, - II V^(1)(0STF + ^>yl/2 ||2^^< ^27/5^-14/5^ (2.86) 

□ 

In the case of TZs^Qi) we proceed as above, namely introduce the auxiliary density 
p(x) = ^[0^'^^(x) + Xn]+'^ and decompose 7^3 (f2i) by using convexity. 

2.14. PROPOSITION. For Z^l'^ < B < Z^ we have 

7^3(fil) = Di,p('\p^^^-p^),^p('\p^^^-p^)) < CB'/'Z'/'. (2.87) 

Proof. By decomposition we have on the one hand the fully semiclassical and 
easier to handle part 

7^4(fil) = D(^«(p^^'^-p),^«(ps^^-p)) < B^ II ^«(P'^"-P) \\l/,< B'/'z'/\ 

(2.88) 

On the other hand there is the more interesting term 

D(^«(p^-p),^«(p^-p)). (2.89) 

For r < we separately calculate 

(f )2D(V'i')e(z, A^), ij?e{z, z- H^,, Xn)) 

< CB^ II ^i')e(z, H^^An) Wl/^ (2.90) 

and 



;^)2D(^!^)[0^ + A;v]f + A^]f ) < CB' II ^«[0^ + A^]f 11^/5 . 

(2.91) 

Whereas ( |2.91| ) can be bounded by 



CB'\\#[Z/\^\Y''\\l/,<CB'/Z\ (2.92) 
( |2.90D can analogously be estimated by [|S[| Lemma 10.7, or [|T| Proposition 4.3. 
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The term 



(2.93) 



is a bit more delicate. We can either use Q Proposition 4.3 or Q Theorem 
4.5.4 (i), which states that for Kh, given in ( p^.40D , with V fulfilhng ( p^.42D and 
|^(x) + r| >e 



e{x,x;Kh,r)~h-'' j 0_(A;(x,p) - r) 



< /i^-'^ Vx G 5(0, 1/2). (2.94) 



Since IAat] < B'^^^Z^^^, we get that |0Stf + Aat] > e/(x)^ in Qi. Hence, we can 
apply (p.94|) to our case, with d = 1, yielding 



\e(z,z;H^^,XN) 



[0^(x) + A;v]f I </(x)-\ 
The term stems from rescaling 5(0, 1) to 5(0,/). So, 



(2.95) 



D{^p!^\p^ - p),^^\p^ - p)) < CB' II 4^^/(x)-i ||^/5< CB'[ZlBf. (2.96) 

□ 

In the case B > Z"^, Z/B is smaller than 1/Z and in the above calculations only 
the separate terms have to be taken into account, which yields analogue estimates 
as above. 



2.5 Analysis in the outer Region ^2 

This region has already been treated by Ivrii in ||l| Section 4. 

Recall first that is the radius of the support of 0^"""^, in the neutral case, 
and of [0^"""^ + otherwise. In order that Theorem p.l| can be applied (f)^"^^ and 
0STF have to fulfill condition ( p:^ ). We know that V(j)^^^{rs) = (p^^^{rs) = 0. 
Hence, we look for a parameter < c < 1, and the concerning radius crs, such 
that the separate quantum mechanical as well as semiclassical parts of TZi-4, 
in {x||x| > crs} do not exceed CB^'^Z^'^ and that 0^^^ fulfills { WM - The 
existence of such a c is a consequence of the behavior of 0^"""^ in the vicinity of 
rs (cf. ( |2.51| )). By means of such a parameter c we decompose the outer region 



^2 into Vl\ U ^2 define a concerning partition of unity, i.e. 

supp-^i = {x| [Z / B]{1 + e) < |x| < cr5}, supp-^g = {x| |x| > crs{l — e)}, 

(2.97) 

with = 1 for r > [Z/ B]{l+e). On i^g? by definition, all terms separately 

are bounded above by CB^^^Z^^^ and on ^2 condition ( p.39| ) is fulfilled for e small 
enough. 

Throughout this section we assume 
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2.15. PROPOSITION. 



(2.98) 



Proof. First we assume B < Z^. Applying Theorem |2.1| with a = 1 and d = 1 
we get for arbitrary but fixed x_|_ (we set z/ = 0) 



IV,:2(M)(#^i(i/xJ)- 



dzdp (2) ,~ , 



< C 



dz/(x)-^/(x) 



(2.99) 



After multiplying with B and integrating over we get 
(2.99) < C— / (ix/(x)-2/(x) < / ^-i^^ < CBZ^/^^^^ji^ ^2.100) 

In the case oi B > Z^ we again have to decompose Ql, since Z/B is smaller than 
1/Z. So for fixed but arbitrary B, fl^ is decomposed with respect to r = 1/Z. For 
r < 1/Z we proceed as in the previous section and estimate each term separately 
and for r > 1/Z we immediately arrive at ( 2.100 ). 
The pure semiclassical part 



B 
2^ 



(ixj_ 



27r 



27r 



can analogously be estimated as in Lemma |2]S . 
Denote 



(2.101) 
□ 



B 



',(2), 



(2.102) 
(2.103) 



2.16. LEMMA. 

i?°(V'f^) < cbIzI. 

Proof. As in Proposition |2.15| we first assume B < Z^. The other case, where 
the terms have to be computed separately, works as in Lemma |2.7|. Applying 



Theorem ^?1] with a = 1 and d = 1 we get for arbitrary but fixed x_l (we set 
iy = 0) 



dzdp , (2) 



2n 



< C 



dzl{^)-^. 
(2.104) 



This implies 

i?0(V'f ^) < C— / rfxZ(x)-^ <CB rdr < CB[rs]'^. (2.105) 
2vr J Jz/B 



□ 



H March 8, 2001 



20 



2.17. PROPOSITION. 

7^2(fi2) < CB^/^Z^/^. (2.106) 

Proof. Note that by the Lemmata PTH] and [^.161 the estimate ( |2.76| ) is proved 
and the assumption of Lemma ^. 10| justified. So by Lemma p.lO| and Lemma |2.16| 



we arrive at ( p.l06| ). □ 

2.18. proposition. 

7^4(^^),7^3(^^^) < cb^/^z^/^ (2.107) 

Proof. Let us start with 

n,{Ql) = D{^p'^\p'^^ - p), ^p^^\p'^^ - p)). (2.108) 
By the HLS inequahty we get 

n,{nl) < CB^\\ij^^\p'^^-p)\\l/, 



< 

The term 



CB^lXj, - Z/|2 II ^f)(0STF + ^)-l/2 ||2^^< C54/5^3/5_ (2.109) 



7e3(n^) = Di^^^\p'^^ - p^), v^f)(pSTF _ p^)) (2.110) 

is a bit more dehcate and we refer to Proposition 4.2 and 4.3 for a proof of 
the estimate (|2.107|) . 

Note: Proposition 4.3 in |^ is proved for region X4 = {^ll^l ^ CqZ/B} 
with possibly a very large parameter Cq. This parameter Co is chosen in a way, 
such that only the lowest Landau band contributes to Ivrii's calculations. Since 
we only treat the lowest Landau band case the assertion of Proposition 4.3 holds 
in our case on the whole region 

Furthermore, we remark that if ( ^.95| ) would be valid on Q^, we could im- 
mediately conclude by the HLS inequality that TZ^i^Ql) < CB'^[rs]^. But since 
the validity of ( p.95| ) cannot be guaranteed on Q\ we have to refer to Ivrii's 
method. □ 

Recall that we have made a partition of unity, Yl'i j=i'^f\'^) — 1- collecting 
all estimations of the Sections 2.4 and 2.5 we have finished the proof of Theorem 
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3 Semiclassical theories approximating ^^^^f 

As we have already argued throughout the introduction, the natural semiclassical 
approximation of -E^^f is given by the DSTF functional 

^DSTF^^j =Y^L I Pm{zf-Z f Vmiz)Pmiz)dz\ +D{p,p). (3.1) 

Here p is a sequence of one-dimensional densities p = (pm(-2))meNo- Iii contrary 
to the usual STF theory the integration over the variables orthogonal to the 
magnetic field is replaced by an expansion in angular momentum eigenf unctions 
in the lowest Landau band. The potentials Vm and D are defined in ( |1.30|) . The 
corresponding energy is given by 



E''^''^{N,Z,B) = mf{S''''^^[p] 

with 



p G P and 5^ /"p„ < iV I , (3.2) 

m ) 



m rn 

Another semiclassical approximation, where the variables, as in the usual STF 
theory, are three dimensional densities, is realized by the MSTF functional 

^""'^^H = ^ / ^^'^'W - / ^xbFV(x) + 5(p,p), (3.4) 
with respective energy 

^MSTF^^^^^^) ^i^f l^MSTFj^] pG^D and y"p<Ar|, (3.5) 

where 

V = {p\ pe £^(R^), j \xpp < oo, D{p, p) < oo}. (3.6) 
First of all, we will show that these two functionals are equivalent. 
3.1. LEMMA. For all peV let 

P(x) = Xm{^±)pm{z), (3.7) 

m 

with pm{z) = J p(x)xm(x_L)(ix_L, ttud denote p = {pm)m- Then one gets 

^MSTF[^]>^MSTF[-]^^DSTF[~]^ (3.8) 
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Proof. By the definition of tlie MSTF functional, it suffices to sliow tliat J > 

For tliis purpose we note tliat for every non- negative function /, on a general 
measure space, one derives from convexity that 

^ '' fdp> ( [ -^F^fdp) . (3.9) 



Hence for every m G N and 2; G M, we have 

-| / p(x)^x™(x^)dx_L > I / p(x)x^(x^)dx_L 1 

\Xm\ Jsuppxm 

Since \xm\ = ^5 we arrive at 

4:7T^ f , , TT^ /■ , ,n , 4:71 



— I p(x)^Xm(x±)c?X_L > / p(x)Xm(x±)rfx_L J • (3.10) 

X,m\ JsuppXm VIXml J J 



352 



p{^fd^ >YJ2 Pmizfdz = ^ / P(x)='f/X. (3.11) 

ra 



□ 



3.2. PROPOSITION. ForallN,Z,B 

^MSTF(^^ Z, 5) = E°STF(;v, Z, 5). (3.12) 

Proof. Lemma p.l| immediately imphes 



E^ST^(iV,Z,i?) = inf{^^ST^[p]|p=^5^p„(^)x™(x^), {p^UeV}. (3.13) 



27r 



□ 



For simphcity we first concentrate on the DSTF functional and then apply our 
results to the MSTF functional. 

3.3. LEMMA. £^^^"'"^[p] is uniformly bounded from below on V. There exists 
even a positive constant a and a C , such that 

>o^{^j pI + Dip, P?j - C (3.14) 

for all p eV. 
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Proof We set p(x) = Zlm Pm(^)l0m(x±)P for an arbitrary {pm)m G 'D. We get 
from ||BBH| Lemma 2 that for every e > there exists a Cs, such that 

j |x|-V < II p II3 +aD{p,py/'. (3.15) 

Hence, this imphes 

J2 / ^™P- <^ 1^/^^ ^5^|0™(x^)|VmWj j +C,D(p,p)i/l (3.16) 

By convexity of x^, for x > 0, and by the equation J2m \'Pm\'^ = we get 

(27r \ 271 

-^5^|0m|Vm <-^5^|0mrpL (3.17) 
m / m 



Using ( 3.171) and integrating over the x^-variable, the inequahty (|3.16|) can be 



C,D(p,p)i/2. (3.18) 



written as 

Consequently the functional £°^"'"^[p] can be estimated from below by (e = 

^°'""[p] > - E / - ^ E / + ^(/^' p) - CsDiP, pf' 

m \ m / 

> inf \kX^ -eX + -CsY] 

X,Y>0 

> a{X^ + Y^) -C, (3.19) 

for a,C appropriately chosen, where we used the notations X = (X]m / PmY^^ 
and y = D(p,p) □ 

3.4. LEMMA. There exists a p^°°' , which minimizes uniquely in T), 

I.e. inf{^DSTF^^]|^ eV]= £°STF[^(oo)]^ 



Proof. Let p*-*-* be a minimizing sequence of £;dstf_ Lemma |3]^ yields that there 
exists a constant C, such that 
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for all m G Nq. By the Banach-Alaoglu theorem there exists a subsequence, still 
denoted as and a p^°°\ with pin G £^(M) Vm G No, such that 

^ weakly in C^{m) \/i G Nq. (3.21) 

Since £^-norms are weakly lower semicontinuous, we derive for all m, 

hminf /(p»)='> /(p(r^)^ (3.22) 

and using Fatou's Lemma we consequently arrive at 

li^L^f E / (Z^™)' ^ E/ (/^™ ^)'- (3.23) 

m m 

Moreover since Vm G £^/^(M) for all m, we conclude by weak convergence 

lim [ dzVUz)pl^{z) ^ [ dzVm{z)p^^\z) (3.24) 
for each m. By (|3.2CI|) and the dominated convergence theorem we have 

/ rf^Kn(^)p^;n^) E/ rf^Kn(^)p^rH^)- (3.25) 
m m 

In order to show 

liminfD(p(^\p») >:D(p(°°),p(°°)), (3.26) 

i— »oo 

we use the fact that for sequences of functions / = (/m(-2))m, 9 = {gm{z))m, 

{f,g)n = Dif,g) (3.27) 

defines a real inner product and consequently a real Hilbert-space T-Cd. Since 
( |OD0 yields || p(*) ||d= V(pW,p»)d < C for all z, we can extract another 
subsequence p^^\ such that 

(/,p»)^(/,p(~))forall/G7^D. (3.28) 

Hence, we conclude 

D(p(~),p(~)) = limi^oo(p^'\p^°°^) < (p(°°\p(°^))^/Miminfi^oo(p^'\p^'^)^/^ 

= ;D(p(°°),p(°°))l/2liminf^^oo^(p^*^P^'0^/^ (3.29) 

and consequently get ( p.26| ). Altogether we have shown 

hminf £°STF[^W] > ^DSTF[^(oo)]_ (3_3Q) 

j— »oo 

The uniqueness follows from the strict convexity of £;dstf_ q 
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3.5. THEOREM. Denote = J pl^\ Then 

(i) for each N < Nc there exists a unique minimizer for S^^'^^ , under 
the restriction I Pm < N, i.e. E^^'^'^{N, Z, B) = S^^^^[p^]. Moreover, p^ 
satisfies J2mIPm = N. 

(a) E^^^^{N, Z, B), as a function of N, is strictly decreasing and strictly convex 
up to Nc, and constant for N > Nc. 



Proof. Let N < Nc- Then the same proof as in Lemma ^.4| shows that there 
exists a p^ eV, with j p^< N and 



Obviously E^^"^^ {N, Z, B), as a function of A^, is non-increasing, and the con- 
vexity of fl^STF i]2iplies the convexity of E^^^^ . Hence, by definition of Nc and 
Lemma |3.4| it is clear that E^^"^^ is strictly decreasing up to A^^^ and constant 
for > Nc. Furthermore we get that j = N ioi N < Nc. (Note that 
Tlim I Pm < ^ would be a contradiction to < Ai'c.) □ 

3.6. PROPOSITION. Let N < Nc. Then for every minimizer p^ there exists 
a parameter p{N), the chemical potential, such that p^ obeys the coupled TF 
equations 

?,K{pi{z)f = [ZVm{z)-Y, I VmAz-z')p^{z')+p{N)]+ V(meNo), (3.32) 

and p{N) fulfills the relation 

^E^^''^{N,Z,B) = p{N). (3.33) 



Proof. The proof works analogously to ||LS[| Theorem 11.10, if the variable perpen- 



dicular to the field is replaced by the angular momentum quantum numbers. □ 



3.7. THEOREM. All statements of Theorem \3.^ and Proposition p. 6| are also 
valid for the MSTF theory, where the minimizing MSTF densities p^(x) and the 
minimizing DSTF densities {Pm{z))m c^re related as 

p''i^) = ^J2xmMpli^)- (3.34) 

m 

The corresponding TF equation reads 

3«:(p^(x))2 = [Z\^^ d^Xm{^A.)Vm,n{z - z')Xn(x'^)p^(x') + p{N)]^. 

n,m 

(3.35) 



H March 8, 2001 



26 



Proof. The existence of a minimizing density p^(x) we get from Theorem ^3 
and Lemma The uniqueness follows from the strict convexity of in p. □ 



Next we try to collect some information about the "critical" particle number 
Nc, which measures the maximal particle number that can be bound to the 
nucleus in the D(M)STF theory. 

3.8. PROPOSITION. > Z. 

Proof. By definition of Nc, we have p{Nc) = 0, so the TF equation reads {p'^" = p) 

2kp„,{z) = (^)]+ Vm G No, (3.36) 

with 

^^^\z) = ZV^{z) VmA^ - z')pn{z'). (3.37) 



We assume < Z. 



BRW| , [RW|| tell us that the potentials Vm{z) and Vn.m{,z — z') behave like l/\z\ 
a.s z —>■ oo. Hence, we get that for each m 



lim \z\[ZVm{z)] = Z, 



(3.38) 



as well as 

lim 1 2; I Vm,n{z - Z')pn{z') 

z~*oo ' / 

- n 

= ^ Jim \z\Vm,n{z - Z')pn{z') = ^ Pn = N^- 

n n 

Since we therefore get 

lim \z\ip''^\z) = Z-N,>0, 

z—*oo 

we can conclude that there exists an e > and a z > 0, such that 

^l^\z) > e\z\ for z > z, 



which by ( p.36|) is a contradiction to pm G ^C^i 



(3.39) 

(3.40) 

(3.41) 
□ 



In the usual STF theory the inequality Nc < Z is a consequence of Newton's 
potential-theory. Since we miss this powerful tool in our DSTF theory we cannot 
expect to get an analogue estimate. But if we use similar methods to those 



BRW 


, Sei 




HSl 



bound for A^^. 



3.9. PROPOSITION. Nr < ^Z . 
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Proof. If we multiply ( |3.36| ) with Pm/Vm and integrate over z, we get 

(3.42) 

Note that by multiplication with pm the []4.-bracket can be dropped, since Pm = 
where (p^^\z) < 0. Clearly J p%_/Vm > 0, so after summing over m we arrive at 

ZNc > / dzdz' . Pn(^OK,m(^ - /)pm(^)- (3.43) 

J ^m(^) 



Moreover, ||HS1|] Lemma 4.1 tells us 



^ + TT^ + T7^ ] Vm.niz - z') > 1, (3.44) 



Vr^iz) VJZ) VrJz') VJZ' 



which we use, together with symmetry, in order to estimate the right side of 



[ \ p.m{z)Vm,n{z - Z')pn{z')dzdz' 

Vm{z) 

m,n ^ ' 

4e 



1111 



Vm[z) Vr,{z) Kn(z') Vn{z' 



Xpm{z)Vm,n{z - z) pn{z')dzdz' > -N^. 



Inserting into ( |3.43 ) finally leads to 



(3.45) 

(3.46) 
□ 



3.10. Remark (The difference between E^^^^ and E^^^). Obviously, the 
magnitude of difference between the D(M)STF and the STF energy is given by 



B 



/,MSTFc^N|3/2 



^STF. x|3/2 



(3.47) 



Due to the singularity of the STF potential, ( p.47|) has to be split into 



B 



/,MSTFc^N|3/2 



|x|<_B-i/2 

+B 



^STF. x|3/2 



|x|>_B-i/2 



^MSTF. ^,3/2 



(3.48) 
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The magnitude of the first term is proportional to 

Z'/^B [ |x|-3/2 = 0{Z''/^B'/^). (3.49) 



The second term of ( |3.48|) could be estimated by 



J|x|>B-i/2 



with 6 = |x_|_|. So we see that the main contribution to ( p.47|) stems from the 
i?^^/^-vicinity of the nucleus, i.e. 

^DSTF _ ^STF ^ o{Z^/^B'/'). (3.51) 

3.1 Some notes about the one-dimensional case 

If we reduce the DSTF functional to the angular momentum channel with m = 0, 
one gets the functional 

£:1DSTF[^] ^ ^ y ^^^(^)3 _z j dzVo{z)p{z) + i y dzdz'Vo^oiz " z') p{z) p{z') , 

(3.52) 

which can be treated analogously to the three dimensional case and Theorem |3.5| 
and Proposition are also valid. Concerning the upper bound of Nc it is not 
necessary to symmetrize over n and m, and in this case ( |3.45| ) reads 



-p{z)Vofi{z — z')p{z')dzdz' 



I n I 1 \ ^ 



V^{z) 

2 / [v^) + V^) ) - -')P(-')C^-C^^' > ^^Nl 

Consequently one gets Nc < 2Z for the maximum particle number that can be 
bound to the nucleus in the one-dimensional theory. 

Moreover, let us regard the absolute minimum E^^^'^^{Z, B) of the functional 
^iDSTF|-^j _ £;^^|^pj^ jf ^f^Q ijigg scahng relations 

V,{z) = B^/'V^{B^/'z), V,,,{z) = B^/'Vl,{B^/'z) (3.53) 

and define 

p{z) = B^'^Z^'^p{B^'''z) (3.54) 

we get 

£zAp\ = B"'Z^I'£lM (3.55) 
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with 



£^^^[p]=K j dzp{zf - j dzV^{z)p{z) + j j dzdzVl^iz - z')p{z)p{z'), (3.56) 



and A = 2B^'^Z^'^. 

Let EIj^{Z, B) be the minimum of the functional 



£i''[p] = K / dzp{zf - Z / dzV^{z)p{z) 



(3.57) 



where the repulsive energy term is omitted. Using the above scaling ( p.54| ) one 
immediately gets ElP{Z,B) = Z^''^B^/^El^{l, 1). So we can formulate the fol- 
lowing theorem: 



3.11. THEOREM. IfZ,B>l are fixed, then 

.iDSTF^^^ 5) < El^{Z, B) + Z{l + 2 \n\BZ^f) 



El''iZ,B)<E' 



(3.58) 



Proof. The lower bound is obvious. 

For the upper bound we use the relation ( |3.55| ) and take the TF-solution of 
£^,, i.e. 

p{z) = -Jv^). (3.59) 

This density is neither in nor in C^, so we define a cut-off density Pr{z) = 
7r^yV^\z)&{R — \z\) and use this as comparison density in ( ^.55| ), which leads to 

^i'iM=^^''(l,l)+ J^(^oi^)f' + \ J PR{z)Vl,{z~z')pR{z')dzdz'. (3.60) 

Since Vqq{z) < min{^, ^/^^J4:} we get by Young's inequality 

I ph) ■^ + 21n(/3) I pI V/5>1. 

(3.61) 

After estimating J p^ and J pj^ we see that the minimum of ( p.61|) as a function 
of P is achieved for f3 = R/\n{R), which implies 



PR{z)Vf^Q{z - z')pR{z)dzdz' < 



Pr{z)V^,{z - z')pR{z')dzdz' < [\n{R) + \n{Rf 



(3.62) 



Next, optimizing the last two term on the right side of ( p.60|) with respect to R 
and multiplying with B^^'^Z^^'^ yields the statement of the theorem. □ 

By aid of this theorem we can also prove that ^^^^"'"^(Z, i?) is the semiclassical 
approximation of Tr£2(]g)[— 9^ — ZVq{z)]^, the sum of all negative eigenvalues of 
-d^.-ZVJz). 
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3.12. COROLLARY. Let B, Z > 1 and B < Z^. Then there exists a constant 
C , such that 

|Tr^2(iR)[-5'-^K)(^)]--^'°^^^(^,5)| < Cmax{Zln[5Z2], 53/4^1/2}. (3.53) 

Proof. This is an immediate consequence of Theorem |3.11| and |Q Theorem 3.19, 
which says that 

[^^cHmi-dl - ZVo{z)]^ - E'J'{Z,B)\ < CB'/'Z'/'. (3.64) 

□ 

We learn from Theorem 3.10 that in a model of a one dimensional semiclassical 
atom, where the electrons are forced to stay in the angular momentum channel 
m = 0, the repulsive interaction energy does not contribute to the leading order 
of the energy E^^^^^{Z, B) for large Z and > 1. 

An analogue effect one obtains for the quantum mechanical interaction energy 
of N particles reduced to the angular momentum m = 0, i.e. 

^ = <i)Q® ...®<pQij{zi,...ZN). (3.65) 

For ip a Slater-determinant or at least for close to the ground state of the cor- 
responding A-particle Hamiltonian iJo? which is the projection onto the angular 
momentum eigenspace with angular momentum m = 0, the interaction energy 
can be bounded from above by (for a precise lower bound see |[HS2|| ) 

^/ p^{z)p^{z')VQfi{z - z)dzdz , (3.66) 

which can be estimated by an analogue method to ( |3.61| ). This leads to 

\ [ p4^)p4^')Vo,oiz - z')dzdz' < CEI'^N"'' [I + In(SAVEo)] , (3.67) 

where we have used that (\E', Hisf"^) < and Eq is the corresponding ground 
state energy (of wave functions of the form ( p.65|) ), which is of the same order as 
^iDSTF long as B < Z^. Relation (|3.671 ) yields that the quantum mechanical 
interaction energy in one dimension is <^ Eq as long as Eq ^ N. 
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